Nonlinearity and optical bistability caused by local 
; field effects in the interaction between a two level 
O atomic system and a metamaterial 

(N : 

\ J L Garcia-Pomar, R Beigang, M Rahm 

I Department of Physics and Research Center OPTIMAS, University of Kaiserslautern, 

■ D-67663 Kaiserslautern, Germany 

, Fraunhofer Institute for Physical Measurement Techniques IPM, D-79110 Freiburg, 

r^i ' Germany 

^ I E-mail: garciapomar@physik.uni-kl.de 

O ' Abstract. We present a theoretical analysis of the nonlinear response of a 

C/3 . metamaterial consisting of an array of split ring resonators coupled to a two level atomic 

I system (2LS) without inversion. We calculated the total complex susceptibility and 
the linear, third and fifth-order nonlinear solution of the Maxwell Bloch equations that 

(— I , describe the coupled system. We analyzed the impact of the nonlinear effects on the 

Oh' dynamics of an active metamaterial/2LS compound in dependence of the magnitude 

of the external electric field when the 2LS is pumped by an external pump source and 

T— I I population inversion is achieved. We observed that nonlinear damping in the active 

' metamaterial/2LS compound significantly reduces the effect of spiking. Finally, we 
investigated the conditions under which the system can exhibit the phenomenon of 

Ps| I optical bistability. We observed that optical bistability can only be obtained when 

\l ' the atomic volume density of the 2LS attains unphysically high values and thus is not 

. expected in the system under investigation. 
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1. Introduction 

In the last decade, the research field of metamaterials has been strongly driven by the 
desire to control or manipulate electromagnetic waves and to realize optical media with 
electromagnetic properties that cannot be found in conventional or natural materials. In 
the realms of the rapidly developing and highly multidisciplinary field of metamaterials, 
several key research directions emerged including negative index materials [H HI [3], 
hyperlensing cloaks and transformation optics |5l E], chiral metamaterials [3, E], 
among others. To date, most known metamaterials are composed of metallic structures 
that are embedded in a dielectric background medium. Since metals suffer from 
significant losses at higher frequencies - as e. g. at infrared or optical frequencies - 
many successful demonstrations of metamaterial applications in the microwave range 
don't immediately translate to other photonic frequency bands. In fact, the search for 
potential applications of the metamaterial technology at optical frequencies is severely 
obstructed by the aforementioned limitations that are set by the intrinsic loss of the 
metallic compositions. In 2007, Boardman et al. [9] proposed to offset the loss in 
metamaterials by incorporating active (gain) media into metamaterials. Although this 
idea is simple at a first sight, it turned out that the coupling mechanisms between 
gain media and metamaterials are complicated and have to be understood in detail to 
accomplish the desired effect. In the last two years, different groups contributed to the 
theoretical description of metamaterials with gain [TOl [11]. For example. Fang et al. 
[12] described a self-consistent approach to describe the interaction between a four-level 
atomic system and a metamaterial consisting of square split ring resonators. Recently, 
Xiao et al. [13] experimentally demonstrated gain in a negative index metamaterial by 
engineering dye molecules as a gain material into a fishnet structure. 

Another intriguing feature of metamaterials at high power levels is the appearance 
of nonlinear effects at high external fields. The possibility to control the effective 
parameters of a metamaterial by exploiting nonlinearities was first suggested in Refs. 
[T4| [15]. In the work described in [H] a nonlinear response of the metamaterial was 
achieved by incorporating a dielectric with intrinsic nonlinearity into the metamaterial 
structure. It was found that the nonlinear effect of the composed system can exceed the 
intrinsic nonlinearity of the dielectric. The reason for this behavior can be explained 
in terms of a strong enhancement of the local fields in the vicinity of the metallic 
structures which causes a strong coupling between the nonlinear dielectric and the 
metamaterial. Consequently, the nonlinear response of the combined system can be 
much larger than the intrinsic nonlinearity of the dielectric. Nonlinear metamaterials 
were successfully demonstrated in the microwave regime by introducing varactor diodes 
as nonlinear elements in a metamaterial consisting of split ring resonators. It was 
shown that the magnetic resonance could be dynamically tuned by varying the power of 
the incident microwaves due to the nonlinear response of the metamaterial. Recently, 
seminal work was contributed regarding the extension of the theoretical description of 
nonlinear metamaterials [16] and the validity of the theoretical framework was verified 
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Figure 1. Schematic model for (a) a two level system with resonant frequency ujab 
and transversal damping "fab coupled to (b) an array of split ring resonators, with 
resonant frequency and damping losbm ^-nd 7s_r_r, respectively, by means of the local 
fields characterized by the coupling factor £. (c) Front view of a unit cell of the system 
where the active medium is surrounding the metamaterial. 



by experimental demonstrations pTl ITSl fT9] . 

Here we present an analytic model model that describes the coupling between a 
nonlinear two-level atomic system (2LS) and a split-ring resonator (SRR) metamaterial 
by means of a density matrix formulation. We show that the coupling between the 
nonlinear medium and the metamaterial is strongly enhanced by local field effects. We 
studied the interaction for two different regimes. In the first regime, we investigated the 
nonlinear effects in a hybrid system consisting of an SRR metamaterial and a passive 
nonlinear material. We derived a full expression for the complex nonlinear susceptibility 
describing the steady state of the coupled system and discussed the contributions of the 
linear and nonlinear susceptibilities of third and fifth order to the total response of the 
hybrid material. We showed that such a coupled system exhibits optical instabilities 
under certain operation conditions. In the second regime, we studied the temporal 
dynamics of the population inversion and the output power of the hybrid system in 
dependence of the incident electric field when the nonlinear medium provides gain and 
is actively pumped by an additional pump source. We found that the spiking in the 
dynamics of the SRR metamaterial coupled to an active gain medium is strongly reduced 
by nonlinear damping effects. 

2. Basic Model 

Metamaterials are usually composed of periodic arrays of metallic basis structures with 
physical dimensions that are significantly smaller than the wavelength of the propagating 
waves. Dependent on the proper choice of the geometry of the metal structures, 
metamaterials can provoke a strong local confinement of electromagnetic fields in the 
medium. For example, the electric field of split ring resonators is concentrated and 
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dramatically enhanced in the gap of the capacitance. In consequence, metamaterials 
can strongly interact with other media that are placed at locations where the local 
electromagnetic fields are sufficiently high. It is obvious that the local field enhancement 
plays a crucial role for the generation of high electric fields to drive the nonlinear 
polarization in nonlinear media. 

In the following, we studied the interaction of a two level atomic system (2LS) with 
a metallic metamaterial consisting of an array of SRRs. We assumed that the coupling 
between both systems was caused by the electromagnetic interaction of the local fields 
that penetrate the SRRs and the atomic two- level medium (Fig. [1]). The metamaterial 
was composed of a single layer of SRRs of copper. The SRRs were described in the 
dipole approximation whereas the 2LS was treated according to a semi-classical model. 
The theory of such a model has been thoroughly studied in Ref. [10]. However, it 
should be noted that the model described in [10\ was inherently limited to linear effects 
in metamaterials since it only took into account a fixed population inversion in the 
2LS that was independent on the magnitude of the incident electromagnetic fields. 
Therefore, the theoretical approach was not suited for the description of nonlinear effects 
in metamaterials and did not account for saturation effects in the gain medium. Below 
we describe a theoretical approach based on a rigorous density matrix description that 
includes nonlinear effects in a coupled system of a metamaterial and a 2LS and allows 
to calculate the complex nonlinear susceptibility of the combined system. 

To determine the magnitude of the linear and nonlinear response of the hybrid 
system we described the local fields of the SRRs and the 2LS as a sum of the external 
field Eext and the local electric fields E^. Thereby, we assumed that the local electric 
fields E^]^^ and E^j^^ at the location of the 2LS and the SRRs were mutually generated 
by the corresponding polarization fields of the SRRs and the 2LS, respectively. In order 
to take the coupling efficiency between both systems into account, we introduced a 
phenomenological coupling constant (local field constant) C . By this means, we could 
describe the local fields as 



where Psrr and Pat are the macroscopic polarization of the SRR and the 2LS 
respectively, given by: 



= E,,, + CP, 

Ef = Ee,, + CP, 



ab 



SRR 



(1) 

(2) 



Psrr — NsRRfisRRPsRR 

Pab = NabfiabPab- 



(3) 
(4) 



Hereby, Nsrr and Nab are the spatial density of SRRs and atoms, Psrr and pab are the 
dipole transition moments of the SRRs and the 2LS for the transition from the ground 
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state I a) to the excited state \b). Applying the rotating wave approximation (RWA), we 
expressed the slowly varying amplitude of the coherence p as 

p = pexpiiut) (5) 

and the electric field as 

Eext = EextCOs{ut), (6) 

In this description the diagonal elements denote the probabilities that the SRR 
or 2LS are in the upper or lower level and the non-diagonal elements determine the 
coherence between the states \a) and \b) [201 [21]. From the density matrix formalism we 
obtained the Maxwell-Bloch equations for the system and could describe the interaction 
of a 2LS with electromagnetic fields at frequencies that are closely tuned to the resonant 
frequency of the atomic transition by the level population difference w and the slowly 
varying off-diagonal matrix element pab'- 

w = -r,, iw + l)+^ [pl^^E^l^ - Pab^Ef^) + r„ (7) 

pab = -i {Aab - ilab) Pab + i^E^^^W. (8) 

Hereby, Aab = co — ujab is the detuning, Tab, 'Jab, are the longitudinal and transversal 
damping of the 2LS, respectively and Tpump is the external pumping. In the above 
equations, the value of the level population difference is —1 if all electrons are in the 
ground state, for equal population in both levels and +1 for total inversion. 

On the other hand, we derive the density matrix of the SRRs from the Liouville 
equation as 

PSRR = -I {^SRR - nSRR) PSRR + l^—El (9) 

where Asrr = u — usrr is the detuning of the external field frequency u from the 
resonant frequency usrr of the SRR and 'Jsrr is the damping constant of the SRR. 
Solving the equation system composed by equations ([7j), ([H]) and ([9]) we obtained the 
solution for the coherences Psrr, Pab and the level population difference w. By inserting 
the solutions of the above system into Eqs. ([3]) and ([4]) we calculated the effective electric 
susceptibility x of the system as a function of the macroscopic polarization fields as: 

Pab{pab,w) PsRRipSRR,w) 
X = ^ab ^ H VSRR ^ l-LUj 

Here, Vsrr and Vat are the volume of the 2LS and the SRR divided by the total volume 
of the unit cell of the complete system, respectively. Note that Eq.( !TO|) determines the 
total susceptibility of the coupled system including all linear and nonlinear terms that 
describe the response to external electromagnetic fields. 



3. Nonlinearity 



We first investigated the stationary nonlinear behavior of the coupled metamaterial- 
2LS compound for the regime where the 2LS was purely passive and not pumped by 
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an additional external pump source Tpu^p = 0. To understand the nonlinear physical 
effects in the system it is illustrative to represent the effective susceptibility as a power 
series expansion of the external electric field by 

X = X^'^+x'^'^\Ee.t\' + x'^'^\E^.t\' + ... (11) 

where we consider only odd orders since the metamaterial structure is centrosymmetric, 
(as the geometry shown in Figll](b)). To obtain the values of the nonlinear susceptibility 
of different order in Eq. ffTTl) . we analyzed the nonlinearity of the level population 
difference w by deriving the power-series expansion with respect to the relative electric 
field strength x = \Eej.t\'^/\Es\'^, where Eg is the analog of the saturation field 
strengthpO]. Under the condition that Asrr = Aab, the saturation field Es is given 
by 

Assuming x to be a small quantity , we performed a Taylor series expansion of w 
retaining only terms up to the second power. For the resonant frequency u = ojsrr = oUab 
the Taylor expansion takes the form: 

1 I (^.2 I ^, ^, \ ^'^^SRrP-SRr/^'^ ~ ISRR ^i^^2/TQ\ 
W = -l+(7^;, + -fablSRR) ; — ^ -^X+OX (13) 

[lab'ySRR + i^fJ'abfJ'SRR/^) NabNsRR) 

From equation (fT3|) we observe that the main variation in the population inversion w 
stems from the local field effects that are represented by the local field factor C. 

Taking into account that the solutions of the slowly varying amplitude of the 
coherence p, from Eqs. (171), (IHl) and ([9]), are given by |10j 



Y ™ h{AsRR-i'ysRR) ) 



" 7 ■ , ^ (Cp,a^,p.sRRiriYN^t.NsRR (^^) 

^ab Hah ^ W {^SRR-ilSRR) 

^SRR + PabNab^PabPSRR/^ r\ 

PSRR = T : (15) 

^SRR — I'lSRR 

where Vtab = f^abEext/^ and flsRR = PsRREext/^ are the Rabi frequencies for the 2LS 
and the SRR, respectively, we observe from Eq. (IT^ that the change of w involves 
two nonlinear effects: the amplitude change of the susceptibility response and the 
modification of the resonance frequency. 

In a first step, we calculated the total nonlinear susceptibility by the set of equations 
([7]) to ffTOj) . As an example, we plotted the imaginary part (Fig. Et^a)) and real part (Fig. 
El^b)) of the total nonlinear susceptibility in dependence of the frequency for x = 10^^ 
and Nab = 3-10^'^ m~^. We also investigated the contribution of the nonlinear response of 
different orders to the overall response by inserting the Taylor expansion of Eq. ( IT3ll into 
the set of equations f lT^ to f|T5|) . As can be seen from Fig. |2l the nonlinear susceptibility 
reveals a resonant behavior and a pronounced minimum at the resonant frequency. A 
comparison between the linear (blue dashed line), third order (green dotted line) and 
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Figure 2. Imaginary part (a) and real part (b) of the nonlinear susceptibility at 
different wavelengths A of the complete system for a high electric external field |i?ea;t|- 
The insets show the differences between the linear (blue dashed line), third order (green 
dotted line) and complete model (red solid line) with values of ujsrr = ^ab ~ 200GHz, 
f^SRR = • m, ^iab = 6.5 • 10-29C . m, -fsRR = 3.4 • 10i3s-i, = 1.5 • 10"s-i, 

Tab = 10i°s"\ NsRR = 5- 102"to-3 and C = 1.7- lO^Om/F taken from the simulations 
and particular values of Nab — 3 • IO^'^to"'^ and x = 10^^ 



total (red solid line) nonlinear response is shown in the insets of Figs. [2]^a) and[2]^b). 
From the calculations, we concluded that the deviations between the complete model 
and the approximated model were maximal at the resonant frequency and therefore 
in the minimum of the imaginary part of the nonlinear susceptibility. Calculating the 
imaginary part of the total nonlinear susceptibility at the resonant frequency we obtained 

ISRR V wiCfiabfJ'SRR/hyNabNsRR- 'JablSRR 

As can be seen from equation (fT6l) . the total nonlinear susceptibility reaches 
a minimum when the second term in the brackets of 0161) approaches to 1. As a 
further property, the imaginary part of the total susceptibility and thus the absorption 
decreases if the local field coupling, represented by the coupling constant £, increases. 
Furthermore, if the product ■jablsRR is kept constant, the imaginary part of the 
susceptibility at the resonant frequency decreases for increasing damping Jsrr of the 
oscillations in the plasmonic structure. 

In Fig. [3], we identified the validity range of the Taylor expansion in dependence of 
the atomic density Nab and the relative electric field strength x by comparing the values 
of the approximated nonlinear susceptibility to the exact solutions of equation system 
([7]) to (eqn:rhoSRR). Since we observed that the deviation between the approximated 
and the complete model were maximal at the resonant frequency of the nonlinear 
susceptibility as discussed above, we restricted ourselves to calculating and comparing 
the nonlinear susceptibility at the resonant frequency by applying Eq. (fT6i) . As a 
criterion for the validity of the Taylor expansion we required that the deviation of the 
imaginary part of the susceptibility that we calculated by inserting the Taylor expansion 
f[T3|) into f|T6|) . from the exact solutions should be less than 0.1%. 
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Figure 3. Validity ranges for tlie nonlinear susceptibility of different orders in function 
of the dipole density Nab of the two level system and the parameter x = lE^^tl'^ /\Es\'^ . 
The parameters are identical to the parameters chosen in Fig. 2, however Nat and x 
vary. 



Applying this criterion we could distinguish four different regimes as depicted in 
Fig. 131 The yellow area represents the range where only the complete model delivers 
accurate results, i. e. we calculated the total nonlinear susceptibility. As expected, 
the higher order terms in the Taylor expansion play a more pronounced role when the 
electric field approaches the saturation limit. Interestingly, the validity ranges are also 
dependent on the volume density of the SRRs that relates to the mutual interaction 
between the SRRs and the interaction of the SRRs with the surrounding nonlinear 
medium. In the green shaded range in Fig. [3] we included terms up to the fifth order 
in the Taylor expansion of the nonlinear susceptibility. As can be seen, the fifth order 
adds only small corrections to the x^^^ response of the composite material. Over a wide 
range of the electric field strength the overall response of the medium can be described 
by a third order nonlinear process as can be seen from the blue shaded region in Fig. 
El It should be noted that the nonlinearity first increases for higher volume densities 
of the SRRs up to a maximum value (see the minimum value of the linear response) 
and then decreases again. This behavior can be seen when the Taylor expansion of the 
level population difference w in Eq. ( IT3l) is inserted into ( IT6l) . For a volume density 
of the atomic dipoles between values of 10^^ and 10^^ we can observe nonlinearities at 
low field strength. This is especially interesting as to the observation of nonlinearities 
in frequency regimes where only low power levels are obtainable as e. g. the terahertz 
frequency range. 
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4. Impact of Nonlinear Effects on the temporal Dynamics of an actively 
pumped metamaterial 

In this section, we analyzed the temporal dynamics of a hybrid material consisting of 
SRRs and a 2-level atomic system that is actively pumped by external pump fields 
at a pump rate Fq. To accurately account for saturation effects in the population 
transfer from the lower to the upper atomic state we described the effective pump rate 
by Tpump = (1 — 'w)/2To. We note that the the effective pump rate Tpump only refers 
to the process of incoherent pumping, i. e. we only pump the populations between the 
levels I a) and \b), however we do not affect the coherence pab between both states. 

Fig. m shows the evolution of the level population difference w and the output power 
of the hybrid system in dependence of time for different external electric field strengths 
Eext and a constant pump rate Fq = 3.5Faf,. For Eext = 10 V/m, we observe strong 
temporal oscillations of the level population difference and the corresponding output 
power. This typical temporal behavior is well known from laser physics and is refereed 
to as spiking. The term spiking describes the strong oscillation peaks in the cavity of 
a continuous wave laser and the emission of strong laser pulses with high electric field 
amplitudes during the initial turn-on-phase. The effect is caused by the development 
of a high population inversion far above the laser threshold during external pumping. 
The absence of appreciable stimulated emission at the beginning of the pump process 
allows the laser gain to rise substantially beyond the threshold gain. As soon as lasing is 
initiated, the population inversion is transferred into the stimulated emission of strong 
laser radiation within very short time, the gain is depleted and the population inversion 
decreases below the laser threshold. As a result, this single spike in the population 
inversion corresponds to a first strong pulse that emerges from the laser. Continuous 
laser pumping allows a periodic revival of population inversion and high laser gain such 
that we obtain a pulse train of laser pulses with enhanced electric field amplitude. Due 
to damping in the gain medium and many loss mechanisms in the laser cavity the electric 
field amplitude of consecutive pulses is damped and after a certain relaxation time the 
laser process reaches a steady state where the population inversion corresponds to the 
threshold inversion and the lasing is continuous at constant output power. In the steady 
state, the emission of laser pulses does not occur and the laser emits a continuous wave. 

As can be seen from Figs. 111(b)- (d), the spiking in the hybrid metamaterial is 
reduced when the external electric field is increased. At an electric field amplitude 
of 10000 V/m the spiking oscillations are almost totally damped out to a single cycle 
oscillation at the very first start of the amplification process. The damping of the spiking 
is caused by an increase of nonlinear damping for increasing external electric fields. The 
nonlinear response in the hybrid metamaterial acts similar to a nonlinear absorber in a 
laser cavity that reduces the spiking. 



10 



(a) 



(b) 



0.5 r 



-0.5 

-It 

I 

15 
= 10 

I 



0.5 1 1.5 2 2.5 3 
xio'° Time(s) xio" 



J[lllllliilMMAMA^AMMA^MA^^^^" 



0.5 1 1.5 2 2.5 3 
Time(s) xio" 




4 

7 3 

i 2 

1 1 



0.5 1 1.5 2 2.5 3 

-10 

>,iq'|i Time(s) xlO 



0.5 1 1.5 2 2.5 3 

-10 

Time(s) xlO 




Figure 4. Temporal dependence of the population inversion w (top) and power of 
the stimulated emission (bottom) for different external electric fields E^^t (a) 10 V/m, 
(b) 100 V/m, (c) 1000 V/m and (d) 10000 V/m. Same parameters as in Fig. 2 and 
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5. Optical bistability 

Optical bistability describes the specific electromagnetic response of a material that 
supports two (meta-)stable states at a given external field stimulus with an electric field 
amplitude E^xt- Optical bistability was first discussed by Hopf, Bowden and Louisell 
[22] . As meaningful applications, optical bistability can be used for controllable optical 
switching and the realization of memory cells. In a recent paper [23], optical bistability 
was observed in a spaser emerging from the introduction of nonlinearities by a saturable 
absorber. In the metamaterial described in this article, we found that optical bistability 
occurs under certain conditions and can be attributed to local field effects that play a 
crucial role in the coupling between the SRRs and the two-level atomic medium. 

For the investigation of optical bistability we considered the Taylor expansion of 
the level population difference w that takes the form 

+ a2w'^ + aiw + ao = (17) 

In the passive case and for Uab = ^jJsrr and zero detuning Aab,sRR = between the 
transition frequency of state \a) and \b) and the SRR resonant frequency usrr we 
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obtained for the coefficients 02, ai and 

a2 = l + 2^ (18) 

-Do 



where ^0 and Bq are described by 



^0 = - lablSRR (21) 
5o = {Cflabf^SRR/hf NabNsRR (22) 

As can be seen from Eqs. (12T|) and ( 122|) . solely depends on the damping of the 
atomic transition 'jat and the SRR resonance 'jab and -Bq involves the dipole moments 
flab and fisRR of the atomic transition and the SRR, the volume densities Nsrr and iVaf, 
of the SRRs and the atoms and especially the local coupling constant C that represents 
the local field effects in the system. 

As a condition for optical bistability it is required that the three roots of Eq. (fTTll 
are real and that at least two of these solutions are stable in the range between -1 and 
+1 in order to provide physically meaningful multiple solutions. We calculated the roots 
of ((|T71)) by numerical solution of the equation system ( IT71) to ( l22l) . 

In Fig. ini we show the value range of the external field Eg^t and the atomic density 
Nab where the condition for optical bistability is fulfilled. Optical bistability can only 
be obtained for high external electric fields and extremely high volume densities of the 
atoms of the gain medium (Fig. [5t^a)). The bistable range constantly broadens for 
increasing electric fields and atomic densities. We note that such high atomic densities 
cannot be achieved in bulk materials and therefore the investigated system is not suitable 
for obtaining optical bistability for the investigated range of values. 

Fig. Mjo) illustrates the level population difference w in dependence of the applied 
external electric field Eext at Nab = 5 ■ 10^^ m~^. As indicated in the figure, we 
can identify two stable branches at a level population difference w = [0, —0.02] that 
approaches zero and thus the population is equally distributed between the upper and 
lower atomic level and ai w = [—0.48, —0.93] where almost all atoms are in the ground 
state. 



6. Conclusions 



We investigated the nonlinearities arising from the coupling of an array of split ring 
resonators (SRR) with a two level atomic medium (2LS) by introducing an analytic 
model that describes the mutual interaction between the SRRs and the 2LS by local 
field coupling. Based on this model, we calculated the total nonlinear susceptibility of 
the system and extracted the linear, third order and fifth order nonlinear contribution 
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(a) (b) 




Eext(V/m) xio= Eext(V/m) xw 

Figure 5. (a) Atomic density Nab versus external electric field Ee^t- The blue zone 
indicates the value range where optical bistability can be observed. Same parameters 
as in Fig. 3. (b) The three solutions of the population inversion for an atomic density 

7V„fc = 5 • 1025 m-3 

by use of a Taylor expansion. We identified the ranges of the nonhnearities of different 
orders in dependence of the atomic densities of the 2LS and the external electric field. 
Furthermore, we investigated the temporal dynamics of the coupled system under the 
presence of external pumping and population inversion in the 2LS. We demonstrated 
that nonlinear damping significantly reduces the impact of spiking during the built- 
up-time of the steady state in the system. Finally, we calculated the ranges where 
optical bistability could occur in the system under investigation. We determined that 
optical bistability could only be achieved under the condition of unphysically high atomic 
densities and therefore is not expected to play a role in the investigated system. 
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